Weak Chaos and Fractional Dynamics in an Optically Driven Colloidal Ring 
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Three colloidal spheres driven around a ring-like optical trap known as an optical vortex have 
been predicted to undergo periodic collective motion due to their hydrodynamic coupling. In fact, 
the quenched disorder in the optically-implemented potential energy landscape drives a transition 
to instability evolving into microscopic weak chaos with fractional dynamics. As a result, the 
relation between the space-time selfsimilarity of the system's collective transport properties and its 
microscopic weak chaos dynamics is revealed. 



Three identical spheres slowly sedimenting through a 
viscous fluid in two or three dimensions generically tum- 
ble chaotically 1]. When the particles are driven steadily 
around a ring, by contrast, their motion is predicted to 
be purely periodic with the reduction in dimension- 
ality and the imposition of periodic boundary conditions 
Hi effectively eliminating the domain of chaotic dynam- 
ics. In this Letter, we demonstrate experimentally that 
chaos can be restored to this system by imposing a small 
amount of quenched periodic disorder. The hydrodynam- 
ically coupled spheres' collective motions, furthermore, 
are characterized by fractional dynamics. 

Our system consists of colloidal silica spheres 1.58 jiia 
in diameter (Bangs Laboratories, Inc. PS04N) dispersed 
in a layer of water 40 /im thick between a glass microscope 
slide and a cover slip. Three of these spheres are confined 
to a horizontal ring 12 /im in diameter by a single-beam 
optical trap known as an optical vortex 5, 6]. An 
optical vortex is formed by bringing a helical mode of 
light to a focus with a high-numerical-aperture lens. 
The helical mode's wavefronts form an ^-fold helix, where 
I is an integer winding number known as the topologi- 
cal charge. The axial ^-fold screw dislocation results in 
perfect destructive interference along the optical axis, so 
that the beam focuses to a ring of light whose radius is 
proportional to 0, 0| • Polarizable particles are drawn 
up intensity gradients to the bright ring, where they are 
trapped. Each photon in a helical beam, moreover, car- 
ries orbital angular momentum that can be transferred 
to trapped objects 0. This creates a constant torque 
that drives the particles around the ring. 

Our samples are mounted on the stage of a Nikon 
TE-2000U inverted optical microscope, whose objective 
lens (100 X NA 1.4 oil immersion Plan-Apo) is used both 
to project an optical vortex and also to create images 
of trapped objects. We imprint helical phase profiles 
onto the wavefronts of a TEMqo beam (Coherent Verdi, 
A = 532 nm) using a phase-only spatial light modula- 
tor (Hamamatsu X6750 PPM) in the holographic optical 
trapping configuration [ijj, |l^ |l3] • 

Figure nja) shows an optical vortex with = 80 whose 
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FIG. 1: Optically driven colloidal ring, (a) Projected inten- 
sity pattern for an optical vortex with ^ = 80. (b) Video 
microscope image of three colloidal silica spheres trapped on 
the optical vortex, (c) Measured intensity variations around 
the optical vortex's circumference. 



image was obtained by placing a front-surface mirror in 
the microscope's focal plane and collecting the reflected 
light with the objective lens. The bright central spot 
is a conventional optical tweezer formed on the optical 
axis by the undiffracted portion of the laser beam. After 
adaptively correcting for residual aberrations in the op- 
tical train 0, , the optical vortex is a nearly uniform 
ring of light with radius Ri = 12 /im. Figure ^b) shows 
three colloidal spheres trapped on the ring at 2.5 W. Un- 
der these conditions, the particles circulate once around 
the ring in T = 1.2 s. We track the spheres' angu- 
lar positions, ^i(t), ^2(^) and ^3(t), by recording the 
video stream at 30 frames per second to a Pioneer 520H- 
S digital video recorder (DVR) with a NEC TL324AII 
monochrome video camera and extracting their instanta- 
neous positions at 20 nm resolution using standard meth- 
ods of digital video microscopy Til. 

The optical vortex's peak intensity, plotted in Fig.^c) 
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FIG. 2: (a) Projection of the phase space reconstructed from 
particle tracks at ^ = 50. (b) Recurrence time distribution 
for ^ = 50. The curve is a guide to the eye suggesting a cycle 
period of 0.6 s (c) Phase space dX I — 80. (d) Recurrence time 
distribution for £ = 80. 

varies by about 20 percent from the mean around its 
circumference. These intensity variations establish an 
effective potential energy landscape for the circulating 
particles resulting from variations in the local orbital an- 
gular momentum flux and from optical gradient forces 
None of these variations is severe enough to trap 
a particle, and speed fluctuations for a single fluctu- 
ating particle suggest that relative mean-square varia- 
tions in the effective potential are smaller than 5 per- 
cent Structure in the potential energy land- 
scape is flxed in both position and time and therefore 
contributes quenched disorder to the particles' dynam- 
ics. 

A circulating sphere entrains flows in the surrounding 
fluid that exert forces on the neighboring spheres. Treat- 
ing the hydro dynamic coupling in the stokeslet approxi- 
mation 19] reveals that the most symmetric state, with 
the particles equally spaced around the ring, is linearly 
unstable against a dynamical state in which two closely 
spaced spheres outrun the third 2]. This periodic dy- 
namical state ought to be observed almost exclusively. 
In practice, however, the experimental system can spend 
much of its time in the nominally unstable state. 

We characterized the system's dynamics by recon- 
structing its phase space from measurements of the inde- 
pendent angular separations, A6i{t) = 02(t) — Oi{t) and 
A02{t) = Os{t) — ^2(^)7 and their associated momenta. 
The results in Fig. El^a) show one continuous trajectory 
obtained over 2.5 hours, and smoothed by box-car aver- 
aging over ten mean circulation periods, T. Smoothing 



suppresses detailed structure in the phase space trajec- 
tory due to diffusion and disorder, and thus provides a 
clearer picture of the system's intrinsic behavior. 

Figure Elja) shows results from an optical vortex with 
£ = 50 and Ri = 8 /im that reveal the predicted peri- 
odic orbit 2]. The parametric trace resembles a torus, 
with a pair of spheres racing around the ring to catch 
up to the remaining isolated sphere to form a transient 
three-particle cluster, whose leading pair then leaves the 
remaining sphere behind. The characteristic cycle time 
of r = 0.6 s substantially exceeds the 0.3 s required for 
an individual particle to circulate once around the ring. 

Increasing the topological charge to £ = 80 increases 
optical vortex's radius and reduces its mean intensity pro- 
portionately. Both the increased interparticle separation 
and the decreased mean circulation rate weaken the inter- 
particle hydrodynamic coupling. The quenched disorder 
is proportional to the local intensity, and thus is com- 
paratively stronger in the larger optical vortex. Under 
these circumstances, the particles' trajectories, plotted in 
Fig. I^c), make frequent and extended excursions away 
from the nominally stable period cycle into the region 
of phase space corresponding to the nominally unstable 
equilateral conflguration. 

We characterized this model system's microscopic dy- 
namics by measuring the distribution of recurrence times, 
r required for a trajectory to revisit regions of phase 
space S6 = 0.5 rad on a side. The results for = 50 
are plotted in Fig. [2fb) and are consistent with a simple 
cycle with a period of 0.7 s. Those for £ = 80, plotted in 
Fig-E^d), display no periodicity, nor are they consistent 
with the exponential distribution that would be expected 
for fully developed chaos. Instead, they follow a power- 
law decay, with an exponent 7 = 1.7 ± 0.1. 

Power-law divergence of trajectories in phase space is 
a deflning characteristic of weak chaos 20]. Its emer- 
gence in periodically kicked overdamped dynamical sys- 
tems has been inferred from measurements on nonlinear 
Alfven waves in plasmas I21L l22l and in low-dimensional 
wall flow of viscous fluids 23]. Unlike these pioneering 
studies, we have direct experimental access to the rel- 
evant microscopic degrees of freedom, and therefore do 
not have to reconstruct the phase space structure from 
measurements of collective properties. 

Quenched periodic disorder is known to induce transi- 
tions to chaos in a wide variety of Hamiltonian and non- 
Hamiltonian systems. The appearance of weak chaos in 
a strongly overdamped system is far less common. Other 
overdamped systems displaying noise-induced chaos Q 
rely on thermal forces to fully explore phase space. The 
optically-driven colloidal ring, by contrast, is athermal. 

Particles circulating around an optical vortex are af- 
fected by three factors: the spheres' hydrodynamic inter- 
actions, the periodic potential landscape, and Brownian 
motion. The importance of each can be inferred by com- 
paring the measured J^] free-particle self-diffusion coeffi- 
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FIG. 3: (a) Evolution of the three- fold bond-orientational 
order parameter, (b) Phase space reconstructed from '^(t). 
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cient 1^0 = 0.19±0.0^um^/s with that of a single sphere 
in an optical vortex Di = 1.0 ±0.2 /im^/s, and with 
that of one of the three spheres in the complete system, 
Ds = 2.0 ± 0.3 /im^/s. The individual sphere's effective 
diffusivity exceeds the free-particle value by a factor of 
five, demonstrating that the landscape's quenched dis- 
order has far more influence than thermal forces. The 
fluctuations, moreover, have a non- Gaussian distribution 
that reflects the intensity variations around the optical 
vortex, rather than thermal forcing. This differs from 
giant amplification of thermal fluctuations that can oc- 
cur when particles become marginally trapped in local 
potential energy wells 11^ • That fluctuations in- 
crease by another factor of two with three particles on 
the same ring demonstrates that hydrodynamic coupling 
is stronger still. Because hydrodynamic coupling alone 
would yield a periodic orbit 2], the weakly chaotic dy- 
namics we observe must be due to the quenched disorder. 
Thermal fluctuations do not contribute appreciably to 
the microscopic dynamics of the observable three-particle 
system, at least in this range of parameters. 

It has long been known that periodic perturbations can 
induce chaos in dissipative systems possessing a stable 
limit cycle j22|. Here, hydrodynamic coupling is respon- 
sible for the limit cycle and the combination of vis- 
cous damping and quenched disorder causes the system 
to sample other parts of phase space. Not any landscape 
will do this, however. As suggested in |22j and proved in 
[2^ , a periodic perturbation can open up a chaotic 
attractor near a stable limit cycle if the perturbation pos- 
sesses a fairly broad frequency spectrum, for instance if 
it consists of sharp kicks. The potential energy land- 
scape inferred from Fig. ^c) has this property, kicking 
the three-particle system three times during each cycle. 

Having direct access to a chaotic system's microscopic 
degrees of freedom presents a rare opportunity to track 
how weak chaos in the microscopic dynamics of an ex- 
perimental system influences macroscopic properties. To 
show this, we combine the three trajectories into the com- 



FIG. 4: Running sum, S{t), of the three- fold bond- 
orientational order parameter, displaying a hierarchy of jumps 
AS spanning a range of durations At. 
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where Oj is measured with respect to a fixed reference 
direction. This function's phase tracks the system's ro- 
tation about the optical axis and is useful for measuring 
the mean cycle period T. Its magnitude, reaches 
unity when the spheres are evenly spaced and drops to 
roughly one third when a pair of spheres is diametrically 
opposite from the third. 

The trace of |'0(t)| in Fig. El^a) is computed for the 
trajectory data in Fig. I2fc) and provides a macroscopic 
overview of the system's microscopic dynamics, akin to 
observing the Brownian motion of a colloidal sphere as 
a probe of the microscopic dynamics of the surrounding 
fluid. As for other effective macroscopic descriptors, 
can be used to reconstruct the underlying microscopic 
phase space. For example. Fig. I2fb) shows a Poincare 
section at delay 20 T, which effectively fills the accessi- 
ble part of the phase space. The periodic state concen- 
trates the trajectory around (1/3, 1/3), while occasional 
excursions to the equilateral state fill out the pattern. 

The order parameter's magnitude also reveals that the 
system switches intermittently between paired and equi- 
lateral states. To analyze the intermittent dynamics we 
use methods originally developed to interpret de nsity 
fluctuation data obtained from tokamak plasmas |29j . 
Given measurements of ilj{t) at N discrete times t/c, the 
running sum 
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minimizes the noise in \ip{t)\ and emphasizes jumps at 
different scales, as shown in Fig. 01 Jumps in S{t) may 
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FIG. 5: Scaling of the magnitude and duration distributions 
for jumps in S{t). 



be identified as fiights in a system with fractional dynam- 
ics. One then can consider the distribution functions, 
Ps{AS) and Pt(At), for the magnitudes, AS, and dura- 
tions At of these fiights. The corresponding results are 
presented in Fig. [SI 

Not only does |'0(t)| display intermittent jumps, but 
its two-state structure appears also to be scale-invariant. 
Both the probability distribution for jump magnitudes 
and jump durations are well described by power laws. 
The exponents a = -2.5 ±0.1 and p = -2.6 ±0.2 for the 
magnitude and duration distributions, respectively, are 
consistent with each other, and both exceed 2. This is 
a signature of fractional kinetics in the system's chaotic 
collective dynamics |2oj . 

Equality of a and /3 arises naturally in any system un- 
dergoing ballistic fiights. For the three colloidal spheres 
in the present model system, these fiights take the form of 
random transitions between the equilateral and periodic 
dynamical states. Furthermore, log-periodic oscillations 
are evident when S{t) is plotted on a logarithmic time 
scale. Such log-periodicityarises for processes character- 
ized by discrete scaling [s^ and can be considered 
as evidence for hierarchical structure in the system's col- 
lective dynamics. The resulting self-similar structure of 
weakly chaotic dynamics is predicted 20] to yield a re- 
lationship between the scaling in the distribution of mi- 
croscopic recurrence times and the distribution of macro- 
scopic fiight durations, P{r) dr = Pt{At) dAt. This leads 
to the prediction 7 = a — 1, which is consistent with our 
experimental results. 

We have demonstrated that imposing a small amount 
of quenched disorder on an optically driven colloidal ring 
can induce a transition from a periodic steady state 
to weakly chaotic dynamics characterized by fractional 
time-space scaling. Because the circulating particles in 
our experiments repeatedly sample the same fixed optical 
intensity pattern, this disorder results in periodic driving 
with strong kicks. Consequently, this provides a model 
system for studying disorder-induced transitions to chaos 
in which both the microscopic and macroscopic degrees 
of freedom are experimentally accessible. 

This work was supported by the National Science 
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